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Abstract. We live near a magnetic star whose cycles of activity are driven by dynamo action 
beneath the surface. In the solar convection zone, rotation couples with plasma motions to build 
highly organized magnetic fields that erupt at the surface and undergo relatively regular cycles 
of polarity reversal. Despite our proximity to the Sun, the nature of its dynamo remains elusive, 
but observations of other solar-type stars show that surface magnetism is a nearly ubiquitous 
feature. In recent time, numerical simulations of convection and dynamo action have taken 
tremendous strides forward. Global-scale organization and cyclic magnetism are being achieved 
by several groups in distinctly different solar and stellar simulations. Here I will talk about 
advances on the numerical front including wreath-building dynamos that may occupy stellar 
convection zones. I will discuss the interplay between the new simulations, various classes of 
mean-field models, and current and upcoming solar and stellar observations. 

1. Introduction 

The 22-year solar activity cycle stands out as one of the most remarkable and enigmatic examples 
of magnetic self-organization in nature. The magnetism we see at the surface as sunspots likely 
originates in the solar convection zone, where turbulent plasma motions couple with rotation and 
magnetic fields to drive strong dynamo action. Magnetism is a ubiquitous feature of solar-type 
stars and there are well known scaling relationships between the amount of surface magnetism 
and stellar properties, such as rotation rate and stellar-type [e.g., Q]]. Other late-type stars 
undergo similar magnetic cycles, with periods ranging from several years to several decades 
0, S|- As in the case of the solar cycle, these must arise through hydromagnetic dynamos 
operating in their convective envelopes. However, the observational landscape is complex, with 
few well-established trends to constrain dynamo models [e.g.,0, Sj- 

Indeed, explaining the origin of globally organized fields and cyclic behavior in the Sun has 
inspired and challenged astrophysical dynamo theory for over a century and continues to do 
so. The tremendous growth of computational resources, coupled with the insights afforded 
by helioseismology about the solar internal structure and differential rotation, has lead to an 
explosion of dynamo modeling efforts. These range from sophisticated 2D mean-field models 
that incorporate data assimilation techniques to fully 3D simulations that can capture non- 
linear dynamics of solar convection and self-consistently establish solar-like differential rotation 
profiles. Both classes of models are being applied to stars other than the Sun, sampling different 
spectral types and exploring how rotation affects stellar convection and dynamo action. 
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Here I will briefly review the state of dynamos for the Sun and sun-like stars that rotate 
somewhat faster, as the Sun did in its youth. These studies are suggesting new modes of global- 
scale dynamo action and are raising exciting questions about the nature of the solar dynamo. 



2. Mean-field models of stellar dynamos 

Simulations of the solar dynamo and solar convection are generally broken into two classes: the 
two dimensional (2D) mean-field models and the three-dimensional (3D) simulations. By going 
to 2D, typically in radius and latitude, the computations can be made much more tractable and 
high resolution simulations can be pursued for very long intervals of time. It is not unusual for 
such models to simulate several tens of solar cycles, and reproducing the entire observational 
record of sunspots is quite feasible [e.g., 0]. 

The trouble lies in the treatment of turbulence and non- linear correlations. Dynamo action 
generally requires correlations between the non-axisymmetric, fluctuating velocity and magnetic 
fields. In mean-field models these non-axisymmetric flows are not directly simulated and instead 
must be captured through some assumed model. This is a very difficult and outstanding 
problem in turbulence theory. Generally, an attempt is made to model the coherent effect of 
the fluctuations in terms of the global-scale (mean) fields, and this is often embodied as an "a- 
effect," though many variations have been explored. Cyclonic convection coupling with rotation 
is often associated with the a-effect and thus a is thought typically to depend on the kinetic 
helicity of the convection, but the a-effect remains very difficult to constrain observationally. 
Comparable difficulties underlie descriptions of the turbulent processes that transport angular 
momentum to self-consistently establish the observed solar differential rotation. This mean 
internal rotation profile, however, can be measured in detail throughout the solar convection 
zone, and the observed profile is directly incorporated in the mean-field models [e.g.,0]. 

The a-effect is most important for the regeneration of poloidal (north-south) magnetic field; 
mean toroidal fields can be generated from a mean poloidal field by the shear of differential 
rotation in what is generally called an "Si-effect." These aSl dynamos and their variants have 
become central to the language used to discuss solar and stellar dynamos. Some live in the 
convection zone alone while others rely on the interface at the base of the convection zone, the 
tachocline, to generate cyclic reversals of global-scale polarity. 

Many recent solar dynamo models have also emphasized the meridional circulation as a 
potentially important factor in promoting cyclic magnetic activity [e.g., @, @-14]. In these 
Flux- Transport models, the equatorward migration of emerging bipolar active regions over the 
course of the solar cycle is attributed to the equatorward advection of toroidal flux in the lower 
convection zone by the mean circulation. Many Flux- Transport models are also Babcock-Leighton 
models whereby the principle source of mean poloidal field generation is the buoyant emergence 
and subsequent dispersal of fibril toroidal flux concentrations, often modeled as a non-local a- 
effect |l5]. Although the physical origin of the Babcock-Leighton mechanism is distinct from 
the turbulent a-effect underlying distributed and interface dynamo models, the fields generated 
are still helical in nature and this helicity ultimately arises from the rotation of the star. 

The review by Charbonneau [TH] is an excellent place to read in more detail about these 
models, while 16f] delves into the details of the turbulence models themselves. Comparisons 
between many of the codes used in different mean-field models have been undertaken by [It}]. 
Mean-field models still do not unambiguously reproduce the solar cycle, despite heroic efforts 
leading up to the current cycle [e.g., @|. To be clear, these models are much further along than 
the 3D models, which are beginning to regularly produce cyclic solutions only in recent time 
[e.g.,[H-[21]. Mean-field models can be run for many solar cycles and are an invaluable tool for 



exploring the extensive parameter space of solar and stellar convection, and several explorations 



have been made of global-scale circulations in solar-like stars [e.g., I22h26||. However, many of 



the important underlying variables are difficult to constrain (e.g., the dependence of a on radius 
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Figure 1. (a) Radial profiles of magnetic diffusivity rj used in various models. Shown is the 
double-step profiles used in the Babcock-Leighton models of [7J and similar to that used in 
271 ] (MF, ). Three profiles from 3D MHD ASH dynamos are shown for comparison; 



these correspond to cases D3, D5 and D10 (grey lines). The molecular diffusivity for a hydrogen 
plasma at solar conditions is also shown, multiplied by 10 7 for display purposes. (6) The Prandtl 
number Pr = v j k (dashed) and magnetic Prandtl number Pm = u/rj (solid) for a hydrogen 
plasma at solar conditions. Simulations use values of order unity. 



and latitude in the Sun; how a varies with stellar properties such as rotation rate and mass; the 
applicability of a-effects to modeling the turbulent induction observed in 3D models; etc.). We 
turn now to a discussion of fully non-linear 3D convection driven dynamos, which are beginning 
to provide the opportunity to better constrain these unknown quantities. 

3. The gap between simulations and the Sun 

Numerical studies of the solar dynamo have a rich history, with the the first 3D 
magnetohydrodynamic (MHD) connective global-scale solar dynamo simulations attaining cyclic 
behavior in the early studies of Gilman [28J. Those Boussinesq simulations were quickly joined 
by fully non-linear global-scale anelastic simulations, which captured the stratified nature of the 
solar convection zone as well [1^, [3(J. Computational resources have grown at a tremendous 
rate, and with them the complexities of the models studied. Modern simulations typically have 
higher resolutions and evolve for longer intervals of time; in one simulation we will examine later 
(case D5), this represents roughly a factor of a million more computation than was possible in 
early studies [28]. This is is in surprisingly good agreement with Moore's law doubling over 
the almost thirty year interval separating these simulations, but it helps clarify the huge gap 
remaining between solar convection and the highest resolution simulations: another century of 
growth might provide the resources to directly simulate solar convection on global-scales. 

Stellar convection spans a vast range of spatial and temporal scales which remain well 
beyond the grasp of direct numerical simulation. Models of stellar dynamos must make various 
tradeoffs, either building up from the diffusive scales or building down from the global-scales. 
These are respectively called local or global simulations; the latter will be our focus here. 
The highest resolution modern 3D simulations, running on massively parallel supercomputers, 
capture roughly 1000 3 total points and typically can evolve for some 10 6 — 10 7 timesteps. In the 
Sun, the smallest scales of motion are set by diffusion and are likely of order 1 mm [e.g., review 
[3l| while the largest scales or motion are comparable to the solar radius (700 Mm), with a total 
spectral range of almost 10 12 in each of three dimensions. Temporal separations are similar, 
with fast granulation on the surface overturning on roughly five minute timescales while the 
deep structure of the Sun evolves over a span of gigayears. 



Clearly, stellar dynamo studies must drastically simplify the physics of the stellar interior. 
Molecular values of the magnetic diffusivity rj for a hydrogen plasma under the conditions of the 
the solar convection zone range from roughly 10 2 -10 6 cm 2 /s as one moves from the tachocline 
to the near surface regions, while the molecular viscosity v is of order 10 cm 2 /s in the solar 



convection zone [e.g., l32l . l33l |. These diffusivities are vastly smaller than the values used in 
either mean-field models or 3D MHD simulations, and we illustrate this in Figure [TJ 

In contrast to the solar values, simulations typically employ values of r\ and v that are of 
order 10 12 cm 2 /s; this large value is more similar to simple estimates of turbulent diffusion 
associated with granulation at the surface where u t ~ VtLt ~ 10 11 cm 2 /s given Vt ~ 1 km/s and 
L t ~ 1 Mm [e.g., l3l| . Shown in Figure [TJa are radial profiles of rj for both mean- field models 
and 3D MHD dynamo simulations with the ASH code. Mean-field models often match to a 
turbulent diffusion consistent with supergranulation at the surface and then taper to a lower 
value at mid-convection zone (here 10 11 cm 2 /s) [e.g., @, 0, [23] - Below the tachocline, t/mf is 
tapered further, sometimes approaching the molecular values. In 3D models the choice is often 
made to scale rj with the background density. Here a scaling of p~ ' 5 is used, though in other 



studies the exponent can scale from to -1 [e.g., [3j[. The relative mixing from diffusion of 
vorticity (by viscosity v), temperature (k) and magnetism (rj) is given by the Prandtl number 
and magnetic Prandtl number, which are shown in Figure Molecular ratios for a hydrogen 
solar plasma are tiny, with Pm ranging from 10 -1 to 10 -5 while Pr is of order 10 -5 . Turbulent 
values are likely of order unity but are not well constrained under solar conditions; simulations 
typically take Pm and Pr to be near unity as resolving large separations in diffusive scales 
requires very high resolutions. 

4. Convection driven dynamos: the Sun 

Despite this daunting separation in parameter space, modern models are making tremendous 
strides in understanding the non-linear couplings between convection and rotation that build 
the solar differential rotation. It is now possible to study relatively high Reynolds number 
convection (fluctuating Re of order a few hundred) in stratified convection zones capturing 
density contrasts exceeding 100 (e.g., more than 5 density scale heights) between the base of the 
convection zone and the near-surface layers. The anelastic spherical harmonic (ASH) code has 
been a very useful tool in global-scale studies 23-3^]; tremendous progress has been made in 



models of photospheric convection as well, but those will not be the focus here [e.g., 1371 ] . 

Simulations of solar convection self-consistently produce solar-like profiles of differential 
rotation, with prograde equators, retrograde polar regions and a monotonic decrease of angular 
velocity with latitude. This profile is achieved partly through the redistribution of angular 
momentum by turbulent Reynolds stresses in the convection 38|-4(J. The differential rotation 



profile is in what is called a "thermal-wind balance" , with an accompanying latitudinal gradient 
of temperature. This is an effect that is well known in the geophysical community, entering 
the vorticity evolution equation as a baroclinic term. In solar convection, thermal-wind balance 



leads to more conical profiles of angular velocity £1 [381 . l40f ] . Simulations suggest that in the 
Sun the accompanying temperature perturbations at the solar surface may be of order 1-10K. 
Perturbations of this size at the base of the convection zone, consistent with the geostrophic 
balances likely achieved in the tachocline, can serve to tilt the contours of constant £1 until very 



good agreement is attained between simulations and helioseismic measurements [411 . [42| . 

Solar dynamo simulations generally produce complex magnetic topologies, with more than 
95% of the magnetic energy in the fluctuating (non-axisymmetric) field components [3]. Mean 
fields are complex, with multipolar structure and transient toroidal ribbons and sheets. Polarity 
reversals of the dipole component occur but they are irregular in time. The presence of an 
overshoot region and a tachocline of rotational shear promotes mean-field generation, producing 
persistent bands of toroidal flux antisymmetric about the equator while strengthening and 



Figure 2. Convection and global-scale flows in case D5 [181 ] . (a) Convective patterns of radial 
velocity v r in global Mollweide projection at 0.95 Rq with upflows light and downflows dark. 
Poles are at top and bottom and thin line denotes stellar surface, (b) Differential rotation shown 
with longitudinally averaged angular velocity fi. The rotation is solar-like, with fast (prograde) 
equator, (c) Meridional circulations with color indicating amplitude and sense of circulation 
(red counter-clockwise; blue clockwise) and mass flux streamlines overlaid. Compared to solar 
simulations, the circulations here are broken into several weaker cells in both radius and latitude. 
(d) Profile of mean temperature fluctuations. This profile, with hot poles and cool mid-latitudes, 
represents the thermal wind balance achieved with the differential rotation, (e) Latitude cuts 

of temperature at fixed radius, sampling top ( ), middle (- -) and bottom ( ) of 

convection zone. At the surface, the temperature contrast in latitude can reach 30K. Profiles 
shown in b—e have been averaged in time over a 225 day interval. 



stabilizing the dipole moment [43-45]. These simulations exhibit notable self-organization 
through the turbulent pumping of magnetic flux into the tachocline, amplification by rotational 
shear, and selective shear-induced dissipation of small-scale magnetic fluctuations. Modern 
simulations build organized magnetic fields and attain cyclic behavior, sometimes with a 
tachocline playing an important role 0, [2(|, and sometimes in the convection zone of the Sun 



alone [21l . 1461 ]. We now turn to the special class of wreath-building convection zone dynamos. 



5. Convection driven dynamos: rapidly rotating suns 

When stars like the Sun are younger they rotate much more rapidly. These stars are observed 
to have strong surface magnetic activity and are thought to have very active dynamos in their 
convection zones. Rotational constraints are stronger in more rapidly rotating systems and this 
can lead to greater correlations as convective structures align with the rotation axis. 

Patterns of convection in a simulation of a young, rapidly rotating sun are shown in Figure Ox 
near the stellar surface. Convection fills the domain and near the equator is strongly north-south 
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Figure 3. Global-scale flows and stellar rotation [47]. (a) Angular velocity shear of differential 
rotation Af2 in latitude near the stellar surface shown as a function of rotation rate relative 
to the solar rotation rate O©. Hydrodynamic cases are shown with diamonds while dynamos 
are labeled and shown with asterisks. AQ grows with more rapid rotation in hydrodynamic 
cases. Cases labeled a and b sample more turbulent states. (6) Kinetic energy of the meridional 
circulations MCKE, normalized to that energy in the simulation at the solar rate. MCKE 
decreases with more rapid rotation; a powerlaw of J7q 9 is shown for reference. Grey dashed 
lines indicate scalings typically used in mean- field models for more rapidly rotating suns (e.g., 
jlsMHol] . and see [5l[ for models that follow the Oq ° 9 scaling). 



aligned. Correlations in these "bannana-cells" transport angular momentum and build the profile 
of differential rotation shown in Figure [2b-; the equator is fast, the poles are slow and the angular 
velocity contrast is larger than in the Sun. The meridional circulations in contrast are weak and 
multi-celled in both radius and latitude (Fig. [2b). Accompanying the angular velocity profile is 
a large latitudinal gradient of temperature, shown Figures [2}i, e. Near the surface, there can be 
30K contrasts between the hotter poles and cooler mid-latitudes. 

The angular velocity contrast of the differential rotation AO and the kinetic energy contained 
in the meridional circulations (MCKE) is shown for many rapidly rotating suns in Figur e CO 
Generally, we find that Af2 grows with more rapid rotation, while MCKE drops strongly [471 ] . 
The decrease of MCKE is a surprise and may hold important implications for flux-transport 
dynamos in the mean- field framework (5ll j. The growth of Af2 is roughly in agreement with 
observations of surface differential rotation in other stars, though substantial disagreement 



remains between different groups of observers [e.g., [52j, [53(. In hydrodynamic cases (triangles) 
this shear continues to grow with faster rotation; in dynamo cases it may begin to saturate 
as Lorentz forces become important and react back on the differential rotation (asterisks). In 
all cases, the growth of Af2 with rotation rate f^o is accompanied by a growing latitudinal 
temperature contrast. The temperature contrast grows from a few K at the solar rate to several 
hundred K at the fastest rotation rates. 

The magnetic fields created in dynamo simulations of rapidly rotating suns are organized on 
global-scales into banded wreath-like structures. These are shown for a dynamo at three times 



the solar rotation rate (case D3) in Figure 54j, [55(. Two such wreaths are visible in the 



equatorial region, spanning the depth of the convection zone and latitudes from roughly ±30°. 
The dominant component of the magnetism is the longitudinal field Ba, and the two wreaths 
have opposite polarities (red, positive; blue, negative). An even more rapidly rotating dynamo 
(case D5) is shown in Figure Hfr. Now the wreaths fill the convection zone and the polar caps. 
These wreaths show significant time- variation and undergo quasi-regular polarity reversals [la ]. 

The time history of these two cases is shown in Figures [SJa, b (case D3 and D5 respectively). 
Here the mean (axisymmetric) longitudinal magnetic field {B^) is shown at mid-convection zone 
over an interval of about 4000 days; both cases have a full history of roughly 20,000 days. 



Figure 4. Magnetic wreaths in stellar convection zones, (a) Persistent wreaths in case D3. Two 
wreaths of opposite polarity (red, positive; blue, negative) form above and below the equator. 
These magnetic structures coexist with the turbulent convection and retain their identity for 
more than 20,000 simulated days. (6) Magnetic wreaths in cyclic case D5. In this simulation 
the wreaths undergo reversals of polarity on roughly a 1500-day timescale. Relic wreaths from 
the previous cycle are visible in the polar caps. This snapshot is at same instant as Figure [2]2. 
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Figure 5. Family of wreath-building dynamo solutions, (a) Time-latitude plot of mean 
(axisymmetric) longitudinal magnetic field at mid-convection zone in persistent case D3 55j ]. 
(b) Cyclic case D5 shown for same span of time [T^]. Three reversals are visible here, occurring 
on roughly 1500 day periods, (c) Primary control parameters for simulations, with magnetic 
diffusivity rj and rotation rate f^o sampling 0.5-15^0. Very approximate dynamo regimes are 
shown. Time dependence emerges at higher magnetic Reynolds number (lower rj) and many 
dynamos undergo repeated reversals of global-scale polarity. Cases with question marks show 
significant time-variation but have not been computed for long enough to definitively establish 
cyclic behavior. At the highest rotation rates the Lorentz force can substantially reduce the 
differential rotation, but dynamo action is still achieved. 



In case D3 we generally find little time variation in the axisymmetric magnetic fields associated 
with the wreaths; small variations are visible on a roughly 500 day timescale, but the two wreaths 
retain their polarities for the entire interval, which is significantly longer than the convective 
overturn time (roughly 10-30 days), the rotation period (9.3 days), or the ohmic diffusion time 
(about 1300 days at mid-convection zone). We refer to the dynamos in this regime as persistent 
wreath-builders. Case D5 (Fig. 0>) is very different. Here global-scale polarity reversals occur 
roughly every 1500 days. Three such reversals are shown here. The ohmic diffusion time in this 
simulation is about 1800 days, while the rotation period is 5.6 days. 

These two simulations are part of a much larger family of wreath-building dynamos, which 
are summarized in Figure Eb. Shown here are 26 simulations at rotation rates ranging from 
0.5 Qq to 15 Q@. Wreath-building dynamos are achieved in most simulations (17), though 
a smaller number do not successfully regenerate their mean poloidal fields (9, indicated with 
crosses). At individual rotation rates (e.g., 30©), further simulations explore the effects of lower 
magnetic diffusivity r] and hence higher magnetic Reynolds numbers. Some of these follow a 
path where the magnetic Prandtl number Pm is fixed at 0.5 (triangles) while others sample 
up to Pm=4 (diamonds). Near the onset of dynamo action the wreaths are similar to those 
found in case D3 and persist for long intervals with little variation in time. At higher magnetic 
Reynolds numbers (lower r/ and higher Qq) we find many simulations that show cyclic reversals 
of global-scale magnetism (as in case D5). It is difficult to determine what sets the cycle period 
in these dynamos: cycles appear to become shorter as r/ decreases, opposite to what might 
be expected if the ohmic time determined the cycle period. The dependence of cycle period 
on Oo is less certain. It is striking that coherent magnetic structures can arise at all in the 
midst of turbulent convection. We find the combination of global-scale spatial organization and 
cyclic behavior fascinating, as these appear to be the first self-consistent 3D convective stellar 
dynamos to achieve such behavior in the bulk of the convection zone rather than relying on a 
stable tachocline of shear. 



6. Where we now stand 

This is an exciting time in solar and stellar dynamo theory. Stellar dynamo models have 
progressed tremendously in the past decade. Mean-field models are reaching a point where 
credible predictions of upcoming solar cycles can be attempted. Meanwhile, several different 
3D simulations using distinctly different codes and assumptions have achieved global-scale 
magnetism and cyclic behavior in simulations of the solar dynamo [e.g., ll9H2ll. |46j, |56J. Global- 
scale organization and cyclic reversals are being found even in simulations without tachoclines. 
A major challenge now is to understand why such cycles occur. Significant progress can be made 
on this problem by translating the results of 3D dynamo models into the language of mean-field 
theory, measuring difficult to constrain quantities such as a and the turbulent electromotive 
force (emf) that builds the mean poloidal fields. First attempts at diagnosing these quantities 



are underway [e.g., [l8|, [55j, [57J but now 3D and mean-field modellers must partner to better 



understand the cyclic dynamo simulations [e.g., l5ll . l58l |. Such efforts will refine the mean- field 



models but will also yield crucial insights into the operation of the dynamos within the highly 
turbulent and time-dependent 3D simulations. 

It is also crucial that further observational constraints be provided for such 3D models; here 
is a brief and biased wish list. In the Sun, some sense of the deep meridional circulations 
would greatly enhance our confidence in modeling results. In particular, it would be useful to 
determine whether these circulations are multi-cellular in latitude and radius, or whether one 
large cell extends from the surface to the tachocline. Simulations generally find that thermal- 
wind balances arise along with the differential rotation. The amplitude of thermal perturbations 
at the surface is likely to be quite small, of order 1-10K, but the detection (or non-detection) of 
such a latitudinal gradient would be very useful. Estimates are beginning to be made based on 



simulations [e.g., 4^, 42l. l59T-[6ll] and attempts have been made to observe this profile in the Sun 
[e.g., [62j, and references therein]. Any estimate of the properties of giant cell convection in the 
Sun would be tremendously useful as well. Direct detection of these structures would of course 
be ideal, but much could be learned from indirect observations as well [e.g., (63[. 

Convective dynamo models in 3D are being applied to other stars, both similar to and different 
from the Sun. For the solar-type stars, observations of the surface differential rotation Af2, its 
variation with stellar mass and rotation rate, and its temporal variation would greatly constrain 
these models. Additionally, it is crucial to know how basic dynamo properties, such as the 
amount of surface magnetism and the cycle period, scale with differential rotation AJ7 rather 
than overall rotation rate fio- The thermal wind balances achieved in more rapidly rotating 
solar-type stars likely lead to latitudinal temperature contrasts of several hundred K; in more 
massive and luminous F-type stars, these contrasts may be as large as several thousand K 
(K. Augustson, private communication). These signatures may be observable. 

The research on wreath-building dynamos in rapidly rotating suns has been done in 
collaboration with Matthew K. Browning, Allan Sacha Brun, Mark S. Miesch, Nicholas J. 
Nelson and Juri Toomre. I owe them a great debt of gratitude. I also thank Ellen Zweibel 
and Cary Forest for inspiring the discussion of rj and plasma transport. Funding for this 
research is provided in part through NSF Astronomy and Astrophysics postdoctoral fellowship 
AST 09-02004. CMSO is supported by NSF grant PHY 08-21899. The simulations were carried 
out with NSF PACI support of NICS, PSC, SDSC, and TACC. 
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